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ABSTRACT 

We study the scenario that a new born strange quark star cools to the Quantum 
ChromoDynamics (QCD) phase transition temperature and converts to a neutron 
star, and we calculate the evolution of temperature and luminosity of the compact 
star. We argue that the conversion energy released can be of the order 10^'^ erg. We 
also propose that a second neutrino burst will be emitted at the completion of this 
phase transition. 
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1 INTRODUCTION 

The Bodmer-Witten proposal that symmetric deconfined 
u,d,s-quark matter may be the absolute g round state of 
matter and fo rms the so-called strange stars jBodmeJIlQTll: 
IWitterJll984^ has aroused much interest, and the proper- 
ties of strange stars have be en vyidely studied ( Alcock 1980; 
ICheng 1998a, : ,Farhl..l984.: Glendenning 1996: Haensel 198S) . 
An important question is whether the observed compact 
stars are neutron stars or strange stars. One possibility to 
distinguish the two is to study their cooling curves, which 
are significa ntly different llWebeilll99gl : In3|2003I : IUsovlll998l 
l200ll.l2002l) . 

In this paper, we study the effects of the Quantum Chro- 
moDynamics (QCD) phase transition on the cooling of a 
compact star and possible signatures of the quark phase. Re- 
gardless of the validity of the Bodmer-Witten proposal, the 
formation of quark-gluon plasma should be favoured in high 
temperature and density l)Shtirva^ll988h : we therefore sug- 
gest that a strange star may be formed just after a supernova 
explo sion, in which both conditions are satisfied ( Bonvcnu tp 
^9£). Because the initial temperature is so high jPetschekl 
19901) Ti ~ 40MeV, the initial compact star is likely to be 



a bare strange star JUsovll200ll) . When it cools down to the 
phase transition temperature Tp, the quark matter becomes 
energetically unstable compared to nuclear matter, and the 
strange star will convert to a neutron star. The conversion 
energy released during this QCD phase transition can be of 
the order 10^^ erg. The temperature drops drastically at the 
completion of the phase transition, which is accompanied by 
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the emission of a second neutrino burst due to the higher 
neutrino emissivity of neutron matter. 

In Section 2, we will first present the equation of state 
(EOS) that we use in the calculation for the quark phase. 
This is then followed by a discussion of the stability proper- 
ties of strange stars in Section 3. Section 4 is devoted to the 
study of the phase transition scenario, which is fixed once 
the EOS's are chosen. We then discuss the cooling processes 
of strange stars in Section 5, and the results of our calcu- 
lation are presented in Section 6. We summarize in Section 
7. 



2 COLD EQUATION OF STATE FROM 
PERTURBATIVE QCD 

Lattice QCD calculation of Tp is highly uncertain at high 
chemical potential, but the latest results llAlltonll2002t) in- 
dicate, though with relatively large uncertainties at high 
chemical potential, that Tp drops from its zero density value 
of 140 MeV to about 50 MeV at 1.5 times nuclear matter 
density po = 0.17 fm~"^ and down to a few MeV for den- 
sity a few times po- Some previous proto-neutron star evo- 
lution calculations indeed show that it is feasible to reach 
the phase transition in supernovae (Benvcrmto 1999). While 
there are still large uncertainties in both high density QCD 
and the proto-neutron star evolution, we believe it is worth- 
while studying the possible consequences of the QCD phase 
transition in supernovae. We assume a constant Tp in the 
star and present results for Tp = 1, 10 MeV for comparison. 
We adopt the simple picture that matter at temperature 
above (below) Tp is in the quark (hadronic) phase. 

To study the properties of quark matter, various EOS's 
have been used. The MIT Bag model is most widely used 
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due to its simp le analytic fo rm llAlcockllQSd) . Here we follow 
Fraga et al. llFragall200il) and use the EOS derived from 
perturbative QCD for cold, dense quark matter up to second 
order in the strong coupling constant as. Qs becomes small 
in the high density limit, with a value of about 0.4 in the 
relevant density regime. It turns out that this EOS is very 
similar t o the MIT B ag model EOS, with an effective Bag 
constant j.Fragal200lf). and we would have obtained basically 
the same results using the latter. None of our results in the 
cooling calculation depends on the validity of perturbative 
QCD in compact star regime. 

All thermodynamic properties can be obtained from the 
thermodynamic potential where /i is the chemical po- 

tential. At zero quark mass limit, the number densities of 
u, d, s quarks are the same, and hence charge neutrality is 
automatically satisfied, without any need of electrons. The 
zero temperature perturbative QCD therm odynamic p oten- 
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iFreedmanlll978l: 



n(/i) = -^^{l-2a, - [G + iV/lna, +/3olnA]Q^},(l) 

where /So = 11 — 2Nf /3, Nf is the number of quark flavors, 
G = Go+Nf{lnNf-0.536), Go = 10.374±0.13, 5, = q^/tt, 
A = A//i, A being the renormalization subtraction point, 
and 
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132 = 2857-325iV|/27, /3i = 2/3i//3gu and P2 = /32/3o/ 8/3?. It 
is be lieved that A lies in the range between 2 and 3 iFraeal 
l20nif) . Both the first and second order terms decrease the 
pressure of the strange quark matter relative to the ideal 
gas. The pressure depends weakly on the strange quark mass 
ms, changing only by less than 5% for nis up to 150 MeV 
JWona 20031 . We wiU therefore use the massless EOS in the 
calculation. 



3 STABILITY OF STRANGE QUARK 
MATTER 

The structure of a static, non-rotating and spherically sym- 
metric strange star can be calculated by solving the Tolman- 
Oppe nheimer-Volkov (T OV) equations together with the 
EOS jGlendennindll99d) . A strange star can be stable even 
at zero temperature if its binding energy is larger than that 
of a neutron star with the same baryonic mass, which is 
indeed the case for A around 2.7 (see T able I), fo r several 
commonly used neutron star EOS's (Webe^El^)- We are 
however interested in the possibility that strange quark mat- 
ter is only stable for T > Tp, and so we choose a A < 2.7, 
so that when the hot strange star cools to low temperature, 
it will convert to a neutron star. For A = 2.473, the max- 
imum gravitational mass is I.5I6M0 with a baryonic mass 
of I.6OM0 and radius 8.54 km. We will use this set of pa- 
rameters in the calculation of the cooling behaviour because 
the maximum mass is close to observational data of compact 
stars. In fact, we have also used other values of A, and the 



cooling behavior is qualitatively similar, as long as the star 
undergoes a phase transition. 



4 PHASE TRANSITION FROM STRANGE 
STARS TO NEUTRON STARS 

It has long been suggested that strange stars can be formed 
from a phase transition of neutron s tars to strange star s 
due to an abrupt increase in density JChenelll99dll998bD . 
However, from the theoretical point of view, formation of 
quark-gluon plasma is favoured whe n both temper ature and 
chemical potential are high enough (|SMr^a3l98i). We pro- 
pose that strange stars are formed in supernovae where both 
the t emperature and density are high, with initial tempera- 
ture iPetscheklll99(]l) Ti ~ 40MeV > Tp. The star then cools 
to Tp and hadronizes into a neutron star containing ordinary 
baryons. This is just the same scenario believed to occur in 
ultra-relativistic heavy-ion collisions |Shurvak 1988). If the 
baryonic mass Mb is conserved during the phase transition, 
the total conversion energy Econv released is: 



= [AfG(SS)-MG(NS)lc 



(3) 



where Mg (SS) and Mg (NS) are the gravitation al masses 
of th e strange star and neutron star respectively llBombaci 
I2OOOI) . Whether a phase transition can occur and how much 
energy is released depend on both the EOS's of quark mat- 
ter and nuclear matter. We choose s everal comm only used 
neutron star EOS's, Ma = 1.4OM0 (IWebeilll99gll . and the 
conversion energy for different A are summarized in Table I. 
For A — 2.473, typically 10^^ erg is released during the con- 
version process, which depends only weakly on the nuclear 
matter EOS. 



5 COOLING PROPERTIES 

The surface of a newborn strange star is so hot that all the 
materials, other than quark matter, are evapor ated leavin g 
the strange star nearly bare without any crust llUsovll200iri . 
Since the thermal conductivity of strange matter is high and 
the density profile of the strange star is very flat, we take 
the uniform temperature and density approximation. The 
strange star cools according to: 

Qf = -L„ (4) 

where Cq is the total heat capacity of all the species in quark 
matter, and Lq is the total luminosity of the star. When 
the temperature drops to Tp, the star undergoes a phase 
transition releasing a conversion energy Econv 

During the phase transition, we assume that the quark 
and neutron matter are distributed uniformly and calculate 
the luminosity of the mixed phase by the weighted average of 
those of the quark matter and the neutron matter llLattimeil 
I1991I: IWondl2003h . When the strange star has converted 
completely to a neutron star, it then follows the standard 
cooling of a neutron star with an initial temperature of Tp. 

The detailed thermal evolution is governed by several 
energy transport equations. We adopt a simple model that 
a neutron star has a uniform temperature core with high 
conductivity and two layers of crust, the inner crust and 
the outer crust, which transport heat not as effectively as 
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the core or quark matter. The typical thickness of the crust 
is ~ 10% of the radius, and we can use the parallel-plane 
approximation to describe the thermal evolution of the in- 
ner crust. The thermal history of the inner crust can be 
described by a heat conduction equation: 

where Ccrust is the specific heat of the inner crust, K is the 
effective thermal conductivity, and ei, is the neutrino emis- 
sivity. As a rule of thumb, the effective surface temperature 
Te and the temperature at the interface of in ner and outer 
crust, Tt,, are related by iGudmundssoniri983h : 

m = 1.288(T,''6/(7.l4)°■"'^ (6) 

where gsi4 is the surface gravity in the unit of 10^* cm s~^, 
Tf,g is the temperature between the inner and outer crusts in 
the unit of 10* K, and Tse is the effective surface temperature 
in the unit of 10*' K. The luminosity at the stellar surface, 
^surface is equal to the heat fiux at the interface of inner 
and outer crusts: 



dT 2 

= -f'surfacc/(47r_R ), 



(7) 



where R is the radius of the star. The boundary condition 
at the interface of the core and inner crust is: 



dT 

'-'core 

Ot 
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(8) 



where Ccorc is the total heat capacity of the core, A^orc is 
the surface area of the core, and L™"^" is the total neutrino 
luminosity of the core. 

5.1 Heat capacity of quark stars 

The total heat capacity is the sum of the heat capaci- 
ties of all species in the star. Without the effect of su- 
perfluidity, the quark matter c an be considere d as a free 
Fermi gas, with a specific heat dlwamotolllOSOl) Cq = 2.5 x 
erg cm K ^ , where p is the baryon density in 
units of po and Tg = T/IO^ K. In the superfluid state, t he 
specific heat is modified as llHorvathlll99il iMaxweHlllQTgTl : 



f 



1.66 + 3.64r 



for T < 1, 



(9) 



where T = T/Tc, keT^ = A/1.76, and A is the energy 
gap in MeV. It has been argued that for quark matter, 
even with unequal quark masses, in the Color-Flavor Locked 
(CFL) phase in which all the three fiavors and colors are 
paired, quark matter is automatically charge neutral and 
no electrons are required ( Rajagop al 2001l"l . However, for 
sufficiently large strange quark mass and the relatively low 
density regime near the star surface, the 2 color-flavor Su- 
perConductor (2SC) phase is preferred. Therefore in a real 
strange star, electrons should be present. The contribution 
of electrons can be parametrized by the electron fraction Ye 
which depends on the model of strange stars. We choose 
Ye = 0.001 as a typical value. The specific heat capac- 
ity of electrons in the strange star is given by fNg 20035 
Ce = 1.7 X 10^° [Yepf'^Tg erg cm"^ K"^ , which is unaf- 
fected by the superfluidity of quark matter. Hence it dom- 
inates the total heat capacity of the strange star when the 
temperature drops below ^ Tc- 



5.2 Luminosity of quark stars 

The total luminosity is the sum of all the energy emis- 
sion mechanisms, including photon and neutrino emission. 
The dominating neutrino emission mechanism is the quark 
URCA process, with emissivity llwamoto 198a) : 



2.2 X 10^^ a, (p) Ye^/^rl erg cm" 



(10) 



and we have chosen as = 0.4 as a constant value throughout 
the quark star. In the superfiuid state, the neutrino emissiv- 
ity is suppressed by a factor of exp( — A/T). 

It has been pointed out that the bare surface of a 
strange star is a powerful source of e^e~ pairs due to the 
strong surface electric fie ld tUsov 19 98). We adopt the e"'"e~ 
pair luminosity given in I^Usmjl2001^ for our calculation. We 
also include the thermal equilibrium and non-equilibrium 
black body radiation in our ca lculation using standard treat- 
ment l| Alcocklll98^ : ll Jso-^ 'ioO 1 ) : the contribution of the lat- 
ter is small compared to other mechanisms at high tempera- 
ture. Once the temperature drops, the cooling is dominated 
by the relatively low power non-equilibrium blackbody ra- 
diation, as long as the star is still in the quark phase. 



5.3 Microphysics of the neutron star 

There are many different models of neutron star cooling. W e 
adopt the one described by Ng l.Ng..200a: .Maxwelll Il979ll . 
The neutri no emission me chanisms are the direct URCA 
processes (iLattimeil Il99ll). the electron-proton Coulomb 
scattering in th e crust iFestal ^G^), and the neutrino 
bremsstrahlung fNe 2003). The surface luminosity will be 
of the blackbody radiation Lbb ~ 47r JJ^ctT*, with the effec- 
tive surface temperature Ts . The blackbody radiation will be 
the dominating cooling mechanism after neutrino emissions 
are switched off. 

For the thermal conductivity of th e inner crust K , 
we use a temperature dependent model jLindblomll2000ll . 
K = 2.8x 10^7^9 erg cm"^ s"^ The choice of K wiU 

not be important after the epoch of thermal relaxation, 
which is of the order 10 — 100 years. The temperature of 
the inner crust and the core will be uniform after that. 
When a strange star is born just after the stellar collapse, 
its te mperature is very high, of the order 10^^ K jPetschekl 
[l99^. A mainly affects the conversion energy, which affects 
the duration of the phase transition. The cooling mecha- 
nisms depend only weakly on A while the cooling curves 
depend weakly on T. We assume T — 40 MeV f or our cal - 
culations. We choose a gap value A = 100 MeV |Npjl20Q^ 
to describe the superfiuidity phase of quark matter, and we 
have checke d that using A = 1 MeV gives qualitatively sim- 
ilar results llWonsll2003 



6 RESULTS 

The observables are the luminosity and the surface temper- 
ature at infinity, L°° and T^ , which are related to the stel- 
lar surface values, L and Ts (Tsuruta 2002): T^ = e'l'^Ts, 
= e^'^'-L, where e'^= = ^1-2A/g/-R is the gravita- 
tional redshift at the stellar surface. The various cooling 
curves (solid lines) for Tp — 1 (10) MeV are shown in left 
(right) panels of Fig. 1. 
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For a large range of parameter values and nuclear mat- 
ter EOS's, we obtain a large energy released, of order 10^^ 
erg. The duration of this energy release depends sensitively 
on Tp ; a higher Tp results in a higher luminosity and shorter 
phase transition duration. It can be as short as seconds for 
Tp = 10 MeV, or as long as hundred thousands of seconds 
for Tp — 1 MeV. The surface temperature of the star drops 
rapidly (top panels of Fig. 1), reaching lO'^ K already within 
the first second {Tp — 10 MeV) to first hundred thousand 
seconds {Tp = 1 MeV). The decrease in surface temperature 
is particularly drastic at the completion of the conversion, 
which is a unique feature of the phase transition not seen in 
either a pure strange star (dotted curves) or a pure neutron 
star (dot-dashed curves). 

The photon luminosity (middle panels in Fig. 1) is ini- 
tially dominated by the e"'"e~ pair emission mechanism and 
is large in the quark phase due to the strong surface elec- 
trostatic field. Since the e'^e~ pairs are not affected by the 
superfluidity gap, the photon luminosity is hardly affected 
by the gap values. The total energy radiated by e+e~ pair 
emission is 4.24 x 10'^°(1.47 x lO'^^) erg for Tp = 10(1) MeV. 
Once the phase transition is completed, the surface field of 
a neutron star is much weaker and this mechanism is turned 
off. The photon luminosity therefore drops by over ten or- 
ders of magnitudes, and this is more drastic for larger Tp, 
because the phase transition occurs earlier and with shorter 
duration. 

A distinct second burst of neutrinos shows up in the 
neutrino luminosity (bottom panels of Fig. 1), which ac- 
companies the completion of the phase transition and arises 
because of the higher neutrino emissivity in neutron matter. 
For Tp = 10 MeV, the neutrino flux rises by over ten orders 
of magnitudes within a small fraction of a second. Note that 
we have not incorporated a detailed transport calculation for 
the neutrinos, which results in a broadening of the neutrino 
bursts we present here. Indeed, for Tp = 10 MeV, the peak 
neutrino intensity of about 10^® ergs~^ lasts only for 10"'' s, 
and these neutrinos wil l be spread out over a diffusion time 
scale of (~ 1 - 10 s) JPetscheklll99ol) . reducing the peak 
intensity by a factor of 10^. The first and second neutrino 
peaks are likely rendered indistinguishable by the relatively 
slow neutrino diffusion out of the dense medium. A much 
more careful treatment of the ne utrino transport is clearly 
needed here jLiebendorfeil 12004V However, if Tp is as low 
as 1 MeV, the two bursts of similar flux can be separated 
by as long as 10^ s, which should be observable by mod- 
ern neutrino observatories. The two neutrino bursts can in 
principle be distinquished also by their energy spectra. The 
flrst burst is emitted near the initial high temperature of the 
newborn strange star, while the second burst is associated 
with the phase transition temperature Tp, and therefore the 
second neutrino burst has a softer energy spectrum. The to- 
tal energy radiated in neutrino is 3.15(1.27) x 10^"' erg for 
Tp = 10(1) MeV. 

This scenario of a second burst of neut rinos can be com- 
pared to two previous similar proposals llBenvenutdll999l : 
lAguileralliooll . In our model, the burst is due to the phase 
transition from a quark star to a neutron star, which has 
a higher neutrino emissivity, whereas in previous propos- 
als, the second burst accompanies the phase transition from 
a neutron star to a quark star. In Benvenuto's theory, the 
phase transition is delayed by a few seconds after the core 



bounce due to the presence of the neutrinos IIBenvenutd 
1999). In Aguilera et al.'s theory iAguilerall20o3) . the burst 
is due to the initial trapping of neutrinos when the temper- 
ature is high and their sudden release when the quark star 
cools. If quark matter is not as stable as nuclear matter at 
low temperature, then there should be yet another phase 
transition back to nuclear matter, which is what we focus 
on, and the "second" neutrino burst we proposed is then the 
"third" neutrino burst. 

If multiple neutrino bursts are observed, as may indeed 
be th e case for the Kamiokande data for SN1987A iHiratal 
Il 98'fl ). whether the compact star changes from the quark 
phase to neutron phase (our model) or the other way around 
can be distinguished observationally in at least two ways. 
First, our model predicts that the cooling is much faster 
before the phase transition, but it will become slower after it. 
Second, the size of the post-phase-transition compact star, 
being a normal neutron star, would be larger in our model. 



7 DISCUSSION AND CONCLUSION 

Based on the lattice QCD phase diagram, we propose that 
the new born compact star in a supernova is a strange star 
and it transforms to a neutron star when it cools down to a 
critical temperature Tp. The conversion energy can be of the 
order 10^"' erg and is adequate to supply energy for GRB's. 
The strange star cools rapidly due to neutrino and e'^e~ 
emission, and its surface temperature drops drastically at 
the completion of the conversion, when a second neutrino 
burst emerges due to the higher neutrino emissivity of neu- 
tron matter. In our models, the phase transition is treated 
in a simplified manner. Hydrodynamic calculation is needed 
for a detailed description of the process starting from a su- 
pernova explosion. Here we discuss semi-quantitatively the 
signatures left by the quark to hadron phase transition, if it 
occurs in supernovae. 

This work is partially supported by a Hong Kong RGC 
Earmarked Grant CUHK4189/97P and a Chinese University 
Direct Grant 2060105. We thank Prof. K. S. Cheng for useful 
discussion. 
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Table 1. Total conversion energy Econv for various A. A neu- 
tron star gravitational mass Mg(NS) = 1.4Mq is assumed, and 
the baryonic masses of strange stars are chosen to equal those of 
the neutron stars. The many-body approximation for HV, HFV, 
A|™P + HFV and G^^|° EOS's are relativistic Hartree, rela- 
tivistic Hartree-Fock, relativistic Brueckner-Hartree-Fock -|- rela- 
tivistic Hartree-Fock and relativistic Hartree respectively (Weber 
1999). 



T71 /~\ O / "\TC ^ 

iijOb(N&) 


Mb/Mq 


A 


71 /f /CON / n ,r 

Mg{SS)/Mq 


£'conv/10'='=' erg 


HV 


1.51 


2.473 


1.44 


+0.72 






2.600 


1.41 


+0.18 






2.880 


1.33 


-1.25 






3.000 


1.29 


-1.97 


HFV 


1.60 


2.473 


1.516 


+2.08 






2.600 


1.478 


+1.40 






2.880 


1.400 









3.000 


1.363 


-0.66 


A|BHF+HFV 


1.62 


2.473 


/ 


/ 




2.600 


1.50 


+1.79 






2.880 


1.41 


+0.18 






3.000 


1.38 


-0.36 


r'K240 
^M78 


1.56 


2.473 


1.48 


+1.43 




2.600 


1.45 


+0.90 






2.880 


1.37 


-0.54 






3.000 


1.33 


-1.25 



This paper has been typeset from a TgX/ X^T^rjX file prepared 
by the author. 
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Figure 1. Cooling curves corresponding to two different mod- 
els. Left (right) panels are for models with small (large) phase 
transition temperature Tp = 1 (10) MeV. The dotted (dashed) 
lines arc the cooling curves of a pure strange (neutron) star with- 
out phase transition. The solid line represents the scenario with 
phase transition. T^, and LJ5° denote surface temperature, 
photon luminosity and neutrino luminosity respectively. Ti = 40 
MeV and A = 100 McV arc assumed. The small insets indicate 
the peaJs luminosities for models with phase trajisition (PT). The 
time axes of the small insets of the upper left, lower left, upper 
right and lower right start at t=155670, 155400, 1 and 1.369 s 
respectively. 



